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Abstract 

We present a calculation of the three-point functions of the 0(A)-symmetric 
sigma model. The calculation is done nonperturbatively by means of a higher- 
order \/N expansion combined with a tachyonic regularization which we pro¬ 
posed in previous publications. We use the results for calculating the standard 
model process f f ^ H ^ WW nonperturbatively in the quartic coupling of 
the scalar sector. 

In previous publications ||I[| we developed a nonperturbative approach for calcu¬ 
lating processes where the quartic self-coupling of the scalar sector of the standard 
model becomes large, and therefore usual perturbation theory becomes unreliable. 
This nonperturbative approach is based on extending the standard Higgs sector to 
an 0(Af)-symmetric sigma model, calculating the scattering amplitudes nonpertur¬ 
batively as a power series in 1/N, and recovering the standard model in the limit 
N = 4. The connection to the physics of electroweak vector bosons is provided by 
the equivalence theorem. 

The idea of expanding in the number of degrees of freedom of the theory under 
consideration instead of the coupling constant is rather old P, ||. This is in principle 
a very attractive idea which attempts to hnd a solution which is valid at strong 
coupling as well as at weak coupling, and which is free of the renormalization scheme 
ambiguities associated with conventional perturbation theory. 
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However, perturbation theory was by far wider used than the 1/N expansion for 
phenomenological purposes. An enormous amount of work in the recent years resulted 
in powerful tools for the computation of higher-loop Feynman graphs. The success of 
perturbation theory is due to the fact that Feynman diagrams can be calculated for 
any theory, and — with increasing difficulty though — in higher loop-orders. 

These are precisely the issues which limited the applicability of the nonperturba- 
tive methods based on the 1/N expansion. First, the precise structure of the 1/N 
coefficients depends on the theory under consideration. This structure can become 
so complicated as to make a direct computation of all the graphs prohibitive even for 
the first term of the 1/N expansion. A classical example of this type is the planar 
QCD in four dimensions [^, where the topological structure of the graphs is known, 
but they could not be actually calculated so far. Second, in most cases of physical 
interest the actual value of N is not that large as to make the leading order a good 
approximation. Such an example is the Higgs physics with which we deal in this let¬ 
ter. Another example is the ordinary QED treated as the A^^e = 1 limit of an Abelian 
theory coupled to N^, species of electrons. When N is not large enough, the leading 
order solution in 1/A^, although simple, is numerically not an approximation of ac¬ 
ceptable accuracy for phenomenological purposes. Thus the inclusion of higher-order 
corrections is mandatory. Higher-order contributions in 1/A^ are however difficult to 
calculate. This is because of technical problems of combinatorial nature, and also 
because of the leak of techniques to calculate certain classes of multiloop graphs. On 
the fundamental side there is the question of treating renormalon-type chains inside 
higher-order 1/N graphs. 

For the case of the 0(A^)-symmetric sigma model, we have shown that the prob¬ 
lems enumerated above can be dealt with, and higher-order calculations in the 1/N 
expansion are feasible iB- The problem of casting all higher-loop graphs into a 
manageable form is solved by using the auxiliary field formalism due to Coleman, 
Jackiw and Politzer 0. The tachyon problem is dealt with by means of a minimal 
tachyonic subtraction [||]. The problem of evaluating the necessary multiloop graphs 
is solved numerically, by adapting a numerical three-loop technique for massive dia¬ 
grams ll, I). 

The two-point functions of the 0{N) sigma model were calculated in refs. |^. 
This leads to a nonperturbative relation between the Higgs mass and width at next- 
to-leading order in 1/N. Since in the scalar sector of the standard model perturbation 
theory has already been extended at two-loop level 0, Si, this allows a strong test 
of the nonperturbative 1/N result at weak coupling. Indeed, at weak coupling there 
is an impressive numerical agreement between two-loop perturbation theory and the 
next-to-leading order 1/N expansion. It is the purpose of this letter to extend these 
results for the three-point functions of this theory. 

As in the case of two-loop functions, we start with the ordinary Lagrangian of 
an 0(A^)-symmetric sigma model. Following ref. 0, we modify this Lagrangian by 
introducing a non-dynamical auxiliary field y: 
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Figure 1: The definition of the subtraeted vertex graphs. The blob on internal lines 
indicates the dressed nonperturbative propagators at LO in 1/N, which result upon 
summation of infinite chains of one-loop bubble diagrams. The wavy line is the a 
field, the solid line is the n field, and the dashed line is the x field. The boxes indicate 
internal and overall ultraviolet subtractions performed at an arbitrary subtraction scale 
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This does not change the physical content of the sigma model, because if one 
eliminates y by using its equation of motion, one recovers the original Lagrangian. 
The advantage of the auxiliary held formalism is that the quartic coupling of the 
<h held disappears, being replaced by trilinear vertices which involve the y held. 
This results into an enormous simplihcation of the possible topologies of multiloop 
diagrams which may appear in higher orders of 1/N. 

The multiloop Feynman diagrams which contribute to the three-point functions 
of this theory at NLO in 1/N are shown in hgure 1. Actually we dehne in this 
picture the subtracted 1/N graphs Ei, Fi, F 2 and F 3 . These are the ultraviolet 
hnite combinations which actually appear in the expressions of observable physical 
quantities upon inclusion of the necessary 1/N counterterms. 

The subtracted three-point graphs of figure 1 are calculated numerically with the 
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same methods as the two-point graphs |^. The main difference is that, because of the 
different kinematic combination of the external momenta, the rotation of the spatial 
component of the loop momentum must be done on a more complicated complex 
path for avoiding the singularities of the integrands. The technical aspects of this 
procedure are discussed in detail in ref. p. 

Once the subtracted 1/N graphs are calculated, they can be used for deriving 
physical amplitudes. Here we will consider the Higgs resonance shape in the scattering 
process ff^H^ WW. This process is for instance relevant for direct searches at 
a possible muon collider. Also it is related to the corrections of enhanced electroweak 
strength to the Higgs production mechanism by gluon fusion at hadron colliders 0 


At NLO in the 1/N expansion, by including the relevant 1/N graphs of all loop- 
orders and the corresponding counterterms derived from the Lagrangian of eq. 1, one 
obtains the following nonperturbative expression for the // ^ if —> WW scattering 
amplitude (the overall factor from the tree level Yukawa coupling not included): 
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Note that this expression is renormalization scheme independent. Here we used 
the notation f}{s) = Ei{s) and 0(s) = Fi{s) -|- F 2 {s) + F^{,s). Similarly, q;(s), /3(s) 
and 7 (s) are the subtracted two-point functions of the model, and 5Z^ and dZ^, are 
the wave function renormalization constants of order 1/Y of the Higgs and Goldstone 
fields, which were defined and calculated in ref. [Q. rn?{s) is the LO self-energy (see 
the notations of ref. [Q), and \fNv = 246 GeV is the vacuum expectation value of 
the Higgs field. 

We plot in figure 2 a set of nonperturbative line shapes of the Higgs resonance 
as it appears in this scattering process. We also plot the location of the peaks of 
these line shapes. For comparison we show also the corresponding location of the 
peaks in perturbation theory, calculated in two-loop order 0,0- Just as in the 
case of the // —;> ff —>• /'/' scattering process []^, the NLO 1/N solution agrees 
well with two-loop perturbation theory at weak coupling. As the coupling increases, 
a saturation effect sets in [|I|, and the width of the resonance increases without the 
position of the peak increasing at the same time. Nonperturbatively, the saturation 
value of the position of the peak is at about 980 GeV. This is not far from the 930 
GeV value found in the fermion scattering process [|^. Of course, the shape of the 
resonance is process dependent, as opposed to the position of the complex pole of the 
Higgs particle, which is universal. 
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To conclude, we calculated nonperturbatively the three-point functions of the 
scalar sector of the standard model by means of a 1/A^ expansion at next-to-leading 
order. Combining the three-point function with the already available two-point func¬ 
tions, we derived the nonperturbative amplitude of the scattering process f f ^ H ^ 
ZZ, WW. Similarly to the already known ff^H^ f'f' scattering, the NLO 1/N 
solution agrees very well with two-loop perturbation theory at weak coupling. At 
strong coupling we confirm the existence of a Higgs mass saturation effect. In this 
process the saturation value is about 980 GeV. This is comparable with the value of 
930 GeV obtained from fermion scattering. 
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Figure 2: The line shape of the Higgs resonance in the scattering process ff^H^ 
ZZ,WW for different values of the guartic coupling. The solid line indicates the 
position of the maxima of the resonances of the nonperturbative 1/N expansion when 
the coupling is increased. The dotted line corresponds to the two-loop perturbative 
scattering amplitude. 
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